O 



Algebras and universal quantum computations 
with higher dimensional systems 

CN ■ 

^ ! Alexander Yu. Vlasov 

! Federal Radiological Center (IRH), 

^— > ■ 197101, Mira Street 8, St .-Petersburg, Russia 

O ■ 

O ■ February 9, 2008 

00 '. 

, Abstract 
> . 

, Here is discussed application of the Weyl pair to construction of uni- 

■ versa! set of quantum gates for high-dimensional quantum system. An 

' application of Lie algebras (Hamiltonians) for construction of universal 

, gates is revisited first. It is shown next, how for quantum computation 

with qubits can be used two-dimensional analog of this Cayley-Weyl ma- 
trix algebras, i.e. Clifford algebras, and discussed well known applications 
to product operator formalism in NMR, Jordan- Wigner construction in 
' fermionic quantum computations. It is introduced universal set of quan- 

^1^, tum gates for higher dimensional system ("qudit"), as some generalization 

I ' > ■ of these models. Finally it is briefly mentioned possible application of such 

' algebraic methods to design of quantum processors (programmable gates 

, arrays) and discussed generalization to quantum computation with con- 

^ ■ tinuous variables. 

: 1 INTRODUCTION 

rS ; 

^ . The models discussed here only recently were introduced in theory of quantum 

" " ' computation H, ^ for description of d-dimensional quantum systems {d > 2), 

but they have rather long history. Initially some special matrices were intro- 
duced by Arthur Cayley more than 100 years ago [Q, and 30 years later similar 
construction was used by Hermann Weyl for building of fundament of quantum 
theory Q . It is so-called Weyl pair of matrices (or elements of infinite dimen- 
sional operator algebra, last case is not considered here in details for simplicity 
and only briefly discussed due to applications to quantum computations with 
continuous variables). The Weyl approach seems very important in theory of 
quantum computation, because it is only way to write an equivalent of Heisen- 
berg relations for finite-dimensional systems used there. In addition algebra has 
important technical applications. It is widely used in quantum computation 
community after reintroducing few years ago with new name "generalized Pauli 
algebra." 
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2 Universal Quantum Gates and Lie Algebras 



"Quantum gate" is special name for arbitrary unitary transformation. "Quan- 
tum networks" |^ may be constructed using many quantum gates and corre- 
sponds to some unitary evolution, "quantum computation." It is reasonable to 
consider some sets of elements that may be used for construction (or approx- 
imation) of any unitary transformation. It is similar with using universal set 
of logical gates (Boolean functions) like AND, NOT in classical networks and 
theory of computations. 

In many models of quantum computations an elementary system is described 
by finite-dimensional Hilbert space, but more difficult continuous case is also 
may be considered. Quantum gates act on some composite system with Hilbert 
space described as tensor product (see Eq. (||) and Eq. (^ below) . For schemes of 
quantum networks each "wire" corresponds to such system and it corresponds 
to exponential growth of dimension of space of states with number of such 
"quantum wires." 

In simplest case the Hilbert space of the elementary quantum system, qubit, 
is two-dimensional. More general case with higher-dimensional systems is also 
used sometimes |^ and considered in present paper. For n qubits Hilbert space 
has dimension 2" and most general quantum evolution may be described by 
2" X 2" unitary matrix. 

The problem of universality is existence of some set or family of unitary 
matrices with possibility of decomposition of any unitary matrix U G U(2") 

U ^U^,U^,---U^^ (1) 

or approximation U w UkiUk2 • • • Ukp with necessary or arbitrary precision, if 
the set is finite k = 1, . . . ,K. 

It should be mentioned, that formally the same gate acting on different 
"wires" corresponds to different matrices, but the additional indexes are omit- 
ted in Eq. (|^) in sake of simplicity. Situation may be more difficult in real 
applications, then system may not be simply presented as tensor productQ but 
it is not considered in present work. 

For quantum computation with n qubits the Hilbert space Ti. has special 
structure due to representation as tensor power of n two-dimensional spaces Ti2 

n 

Due to such structure any unitary transformation of some subsystem with 
m < n qubits, Usuh G U(2™) can be expressed as transformation that acts 
only on m elements (indexes) in tensor product with n terms Eq. (^ used for 
representation of the Hilbert space for n qubits. Such transformations are called 
m-qubit gates. 

^For example in theory of decoherence free subspaces, it is more known, but rather specific 
example. 
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The constructions may be simply extended to higher-dimensional quantum 
systems. Let d > 2 is dimension of Hilbert space TLd for one "qudit" , then n 
systems may be described by d^-dimensional space 

■H^nd®---®nd, (3) 



and unitary group U(c?"). Non-binary quantum m-gates correspond to elements 
of U(d'"), m < 71. 

An important result of quantum theory of computation is existence of uni- 
versal set with two-qubit quantum gates [0, ^, ^ ■ 

There is so-called infinitesimal approach to construction of universal set of 
quantum gates. It was discussed in original papers |0, H, |o| and described 
here only briefly. From physical point of view it can be considered as using 
Hamiltonian approach with Hermitian operators Hf^ (Hamiltonian) , instead of 
S-matrix-like approach with unitary gates (see also Eq. (||) below). 

The infinitesimal approach uses Lie algebra u(iV) of Lie group U(Af) of uni- 
tary matrices!^ together with finite set of elements of this algebra repre- 
sented as some anti- Hermitian matrices = ^^I- It is possible also to write 
iAk = Hk for some Hermitian matrices Hk ~ Hj^. Quantum gate may be 
expressed as 

- cM-^HkT) (4) 

with some real parameter t. Now instead of universal set of elements of Lie 
group V{N) used in Eq. (|l|) it is possible to work with elements Ak — —iHk of 
Lie algebra vi.{N). 

Lemma If the elements Ak generate full Lie algebra u(N) by commutators, 
then it is possible to use set of gates described by Eq. as universal set of 
gates 

Subtleties of transition from Lie algebra to Lie group of gates based on 
Eq. i s not discussed with more detail and may be found in original works 
0, |[ |l0| |2^ , 01 . Description of specific mathematical problems may be found 
also in Ref. 

In physical applications Hk = iAk correspond to Hamiltonian and real pa- 
rameter T, may be infinitesimal [Q, irrational Q or it is possible to consider 
Eq. (m as one-parametric family, where t is time llG] . 



3 Clifford Algebras 

If there is some associative algebra A, then it is possible to introduce structure 
of Lie algebra by using commutators [A, B] = AB — BA. Due to relation 
with universal set of gates discussed above, it is reasonable to look for algebras 
with "simple commutation law" with purpose to have appropriate tools for 
construction and description of universal quantum gates. 

■^Here N = 2"- for quantum computation with n qubits and N = d" for higher-dimensional 
case. 
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As an example of such a commutation law it is possible to mention 
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AB = CBA (yl,BeA C £ C) (5) 

It is commutative algebra for C = 1 ^^nd anticommutative one for C = — 1. It is 
useful also to consider more general case, then C, is arbitrary complex number and 
not all elements of algebra satisfy equation Eq. , but only dijf event generators 
of the algebra. 

Well known examples for C = — 1 are algebra of all 2 x 2 complex matrices 
with generators are Pauli matrices^ 

? J), '^y=(° ~o), '^z=(j (6) 
and algebra of all 4 x 4 complex matrices with generators are four Dirac matrices 

(7) 

Both algebras are Clifford algebras ^ defined by (anti)commutation law 

e^Gj + ejGfe = 2(5j7c, (8) 

where are generators of the algebra. 

Main property of Clifford algebras used in mathematical and physical appli- 
cation is algebraic decomposition of square root 

(^^Cfeefe)' = ^^4, CfeSC. (9) 

Really Eq. (^) often used for definition of Clifford algebra instead of Eq. (^) . For 
example, the property is used for representation of Dirac equation as square root 
of Klein-Gordon equation. On the other hand, Eq. (^) seems not related directly 
with applications of Clifford algebras in quantum computations discussed in 
present paper. 

The Clifford algebras with even number of generators 2n are especially con- 
venient for description of quantum computation with n qubits due to so-called 
product operator formalism [p^ [l^ , [l9| and some other applications [po| [2^ , ^ . 
It is enough to look on standard presentation for generators of such algebra 

[01 

e2k = ii_ 

e2k+l = i 1 ■ • ■ ® 1 (g) (Ty (8) <Tz (g) • ■ • ® (T^ , (10) 

-1 



•^Formally only two Pauli matrices arc enough. Here is some subtlety, because 3D real Clif- 
ford algebra uses all three Pauli matrices as generators and it is not some pure mathematical 



trick, because may be used in description of real NMR experiments [ 



jome Dui 
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where fc = 0, . . . , n — 1 and 1 is unit matrix. 

There are 2^" different products of 2n generators Eq. (10) and it cor- 



responds to presentation of the CUfford algebra via 2" x 2" complex matrices. 
Such representation is useful both for quantum n-gates as elements of Lie group 
U(2") and Hamiltonians as elements of Lie algebra u(2") of that group. 

Jordan- Wigner representation of fermionic creation and annihilation opera- 
tors via ^ ^ 

"I = 2^^2k + ie2fe+i), Ok = ■7^{fi2k - je2fe+i) (11) 
with standard properties 

{afc,aj} = {a[,a]} = 0, {a^,aj} = 4i (12) 



provides useful links with fermionic quantum computation j^, 20, ^| 

It should be mentioned, that an approach with creation and annihilation op- 
erators was used for description of universal quantum simulators and computers 
already in first works of Feynman [|l5[ |l^ , but the question again has some sub- 
tleties, because for description of quantum computation it is also possible to use 
other operators produced via Eq. ( ll|) from slightly changed with 1 instead of 
any cr^ in tensor product Eq. (|l0|). Sometime they are called now parafermionic 
operators, but the special term may be a bit misjudged, because such operators 
even simpler and more natural in formal theory of quantum computation, than 
"usual" fermionic operators. Formally, both kind of operators may be used 
equally in most constructions with creation and annihilation operators used in 
Refs. ll|,|l|. 

Let us return to theory of universal quantum gates. It can be shown, p2t 
that the commutators of elements Eq. ( p^ ) generate only (2n^ -\- n)-dimensional 
space and so set of gates based on these elements via Eq. produces only 
(2n^ + n)D subgroup in 2^"D group of all unitary gates j22|. It is interesting 
counter-example of widespread belief about "omnipresent universality." Any- 
way, appending of only one element of third order like ieoeie2 produces uni- 
versal set of gates [^2| . For fermionic quantum computation sometime is used 
representation with even number of operators Eq. ( pi] ) and so the result is in 
agreement with necessity of elements of second and fourth order for universality 
in such a case . 

The representation Eq. ( [lo| ) is good for general algebraic approach or fermions, 
but in usual theory of quantum computation more useful elements have minimal 
amount of non-unit terms in tensor product. Formally Eq. (^|) produces univer- 
sal set of fc-gates with fc = 1, . . . , n, but it is very simple to construct universal 



*In Ref. |21 1 were considered some points relevant to present consideration, but terms "lo- 
cal" and "spinless fermion" seem not very appropriate. For example, fermionic commutation 
relations Eq. have representation with two-matrices and it corresponds to two-component 
wave vectors of particle with spin-1/2, so spin may not be "omitted." It is also not a model 
of spatially localized particles. If coordinate part of wave function is taken into account, only 
particles with definite momentum satisfy simpler expression like Eq. (h3) and it corresponds 
to maximal "spreading," exp{ikx). 
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set of two-gates. It is enough to consider elements of second order efcefc+i, be- 
cause any such element is one- or two-gate. So universal set of 2n -I- 1 two-gates 
may be based on following elements of Lie algebra u(2") (Hamiltonians) 

Go, Gfc^fc+i = iefcGfc+i (fc = 0, . . . , n - 2), eoi2 = jeoeie2. (13) 

It may be checked, the last element of third order corresponds to one-gate. 



4 Generalization with Weyl Pair 



The construction with Clifford algebra described above may not be applied 
directly for non-binary quantum circuits, when each system described by d- 
dimensional Hilbert space, d > 2. Hilbert space of composite system with n 
qudits is ^"-dimensional and gates belong to U(d"). But even in such a case 
the general principles are very close |^ . 

First, let us consider instead of Pauli matrices Weyl pair ||, Q described 
by Eq. (|) with C = exp{2TTi/d): 



UV ^ exp{27ri/d)VU, VV'^ = C/C/^ 
where U and V are d x d complex matrices [|j 



(14) 
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(15) 



Lately it is called sometime "generalized Pauli group" |l) . 
It is possible to consider three elements 



V, 



with properties 



(16) 



(17) 



Let us introduce set of complex d" x d" matrices similar with construction 
Eq. iM) above: 



= l(g)---(8)l® Tx(8)Tz 



t2fe+i = l®---(Kil®-ry(KiTz 

-k-l 



(18) 



where fc = 0, . . . , n — 1. 
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The 2n elements itself meet to analogue of commutation law Eq. (Q): 

tjtfe = Ctfctj (j < k), {tkf = 1. (19) 

As it was suggested earlier, the rule makes test of universality simpler and the 
elements are really useful for construction of universal quantum gates for higher 
dimensional case (d > 2) 

It is interesting, an analogue of Eq. is also valid ||2^ 

So the elements may be used for expression of algebraic d-th root, but it is not 
relevant to present work. 

The elements Eq. (^8|) are not Hermitian, but for construction of universal 
gates via Eq. (||) it is possible to use |^ 

t+ = z(tfe+4), t- = (t,-t^). (21) 

Commutators of the elements generate whole Lie algebra and so exponents 
Eq. of elements Eq. (|2^) may be used as universal set of quantum gates 

Formally initial expressions Eq. (|18|) represent /c-gates {k ~ 1, . . . ,n), but it 
is again possible to use pairwise products of the elements 



Cfe,fc+1 



= tktk+i {k = Q,...,n-2) (22) 



to construct universal set of two-gates. It is also necessary to use generators 
Eq. with Hermitian operators expressed as sums 



tfc.fc+i - «(tfc.fc+i + tLfc+i)> tfe,fc+i - i'tkM+i - tLfe+i)- (23) 

It should be emphasized, that for higher-dimensional quantum systems with 
d > 2 it is enough only 2n elements Eq. (18) and An elements Eq. ( ^ ) of first 



order and it is not necessary to add some generator of fourth or third order 
for d — 2 (qubits and fermionic quantum computation). Such point may be 
essential, for example elements of high order have difficulties of implementation 
for fermionic quantum computation pl| and even in more general case such 
elements may have specific properties |22| ] . 

If to recall that without such additional element quantum gates generate 
only state of space with quadratic dimension ("complexity") instead of expo- 
nential j2^, |2^, then it is possible to suggest, that for d ~ 2 the question about 
additional element of third order may be essential from point of view of uni- 
versality, correspondence between classical and quantum complexity, etc.. So 
absence of such subtleties with one additional element for higher-dimensional 
system d > 2 maybe not so formal, as it could be suggested at first sight. 

It was mentioned earlier relations of different standard techniques with the- 
ory of quantum computation: unitary gates (S-matrix), Hamiltonian, occupa- 
tion numbers (annihilation and creation operators) . The approach discussed in 
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present section has relation with yet another standard model, Weyl quantiza- 
tion @. The Weyl pair U, V Eq. (|l5|) has many other applications in theory of 
quantum computations |^ , but specific new terminology used in some papers 
sometimes hides it. 

The Eq. ( p^ is Weyl representation of Heisenberg commutation relation. 
The Weyl approach may be even more significant in area or research under 
consideration, because Heisenberg relation itself may not be satisfied in finite- 
dimensional Hilbert spaces used in many applications of quantum computations. 



5 Universal Quantum Processors 

Formally, the quantum network approach used above (and in many other papers) 
resembles classical situation, then some chosen set of logical gates ( "chips" ) are 
fastening on some board to build networks with necessary functions, on the other 
hand, usual electronic computers have fixed structure and different functioning 
is ensured by programmability. 

An analogue of such classical computer is programmable quantum gate ar- 
rays |2^, ^ or, simply, quantum processors |2|, |o|. Here are exists two 
different designs "stochastic" H, ||l and "deterministic" |||, |7[ |2|, |§. 
Formally "deterministic" and "stochastic" processors initially have very different 
design and only recently was found possibility of some confluence in continuous 
limit ||]. 

Anyway only first kind of processors will be discussed in present paper, so 
word "deterministic" is omitted further for simplicity. Such quantum processor 
uses some version of "controlled" or "conditional" quantum dynamics In 
such approach some quantum "wire(s)" are used for encoding of operations on 
second set of "wires." Quantum processor is special network with quantum sys- 
tem represented as tensor product of two subsystem: "program bus" and "data 
bus," but there is some specific difficulties, because different states of "program 
bus" should be orthogonal |^ and so in some meaning all valid programs are 
"pseudo-classical," i.e., a superposition may not be used and would produce 
"malfunction" of quantum processor (entanglement of data and program, more 



details may be found in original works ||25|, |2^, |2^, |2^, |2S|, |30|) 

Let us consider a particular design of quantum processors |g7|; . There are 
three buses: quantum data bus, there any superposition are possible, intermedi- 
ate bus used for choice of particular quantum gate applied to data bus on current 
state. Intermediate bus is also some quantum system, but only fixed number of 
orthogonal basic states may be used to produce proper work of quantum pro- 
cessor. And, finally, there is (pseudo-) classical bus necessary to programming 
of changes in intermediate bus for performing necessary sequence of quantum 
gates on quantum data bus. 

In simplest case the pseudo-classical bus is simply "cyclic ROM," but se- 
quences of universal gates used for implementation of some quantum algorithms 
may be too long for realistic implementations using the ROM based approach. 
So it is better, if sequence of gates may be described by some simple law and 
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processor may use some program of minimal size, i.e. it should be set of in- 
structions like 

repeat { repeat 10 times; repeat Uj 20 times } 3 times, 

but these issues related with application of reversible classical programming are 
out of scope of present work. 

It is clear, that theory of (deterministic) universal quantum processors re- 
sembles theory of universal quantum networks, but sequence of quantum gates 
is not implemented statically, and each next operation is dynamically chosen on 
each step of evolution, accordingly with program. 

At present time a statement, that some particular set of gates is universal, 
has limited significance, because property of universality is already proved for 
very general sets of quantum gates. Discussed algebraic approach is important 
due to constructivity, i.e., there are described concrete and not very difhcult 
algorithms H^, |^ for decomposition or approximation of arbitrary unitary op- 
erators. 

It may be not very actual for some standard quantum algorithms, when con- 
struction of quantum networks is already developed and optimized during many 
years in works of different researchers, but universality of quantum processors 
suggests existence of algorithm for decomposition of arbitrary unitary matrix 
with given precision. Maybe realistic applications are related with mixed ap- 
proach, when the universal algorithm of decomposition is used for fc-gates with 
not very big k (i.e., has appropriate order) and the gates are used as building 
blocks for more difficult quantum algorithms. 



6 Continuous Case 



The Eq. (|14D is Weyl representation of Heisenberg commutation relations and 
valid both in finite and continuous (infinite) case For d ^ oo two families 
of operators with property Eq. (^ may be expressed as 

Ua. Hq) ^ ^J{q + a), Vf,: ^{q) ^ e^'«V(g), (24) 

where a, 6 are two real parameters and C = exp{iab) in Eq. (^. It is possible 
also to write the operators in exponential form 

Ua = e^'^P, Vb = (25) 

where p, q are operators of coordinate and momentum j^. 

All necessary construction used for finite d above may be generalized on 
continuous case |24 , but such universal gates for continuous case have "bi- 
exponential" structure, like 

g{Ua-ul)T _ exp([exp(iap) — exp(— iap)]T) , 
gi{Ua+ul)T _ exp(i[exp(jap) -I- exp(— iap)]T) . 
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Such approach not only produces universal set of gates for one quantum 
continuous variable like in some other works p^ , but for arbitrary amount of 
such variables. It is only necessary instead of tk Eq. in all expressions uses 
their continuous analogues described in Ref. All constructions are based 

on straightforward change of U, V Eq. (|l5|) to Ua, Vt Eq. (p^). Here infinite- 
dimensional Hilbert space is represented as space of complex function with n 
real variables. 

Most difficult part of proof used in Ref. Q related with consideration of 
different nonstandard cases with {^'^ — 1, but it becomes trivial in infinite- 
dimensional case, because it is possible to choose ^ 1 for any integer k. 
Here is again possible to construct set of universal two-gates for arbitrary finite 
number of continuous quantum variables. 

Physical meaning may depend on choice of particular conjugated variables 
p, q. For example, if q is usual coordinate, expressions 

g(v.-yj)r ^ exp(-isin(&a;)T), e*(^''+n^)^ = exp(i cos(&x)r) , 

look more realistic as Hamiltonian for particle in periodic potentials, but more 
difficult analogues of operators Eq. ( p^ ) maybe look more strange. On the other 
hand, the continuous generalizations of are based on some non-canonical form 
of Weyl commutation relations |Q and so may be relevant to description of real 
systems via Weyl quantization (and analogues of Wigncr functions). 
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